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ABSTRACT 

Time evolution of tunneling phenomena proceeding in thermal medium is stud- 
ied using a standard model of environment interaction. A semiclassical probability 
formula for the particle motion in a metastable state of one dimensional system put 
in thermal medium is combined with the formula of quantum penetration factor 
through a potential barrier, to derive the tunneling rate in medium. Effect of en- 
vironment, its influence on time evolution in particular, is clarified in a real-time 
formalism. A nonlinear resonance effect is shown to enhance the tunneling rate at 
finite times of order 2/r7, with 77 the friction coefficient. In the linear approximation 
this effect has relevance to the parametric resonance. This effect enhances the pos- 
sibility of early termination of the cosmological phase transition much prior to the 
typical Hubble time. 
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I Introduction 



Tunneling phenomena, when they occur in isolation, are genuinely a quantum 
effect. But when they occur in some surrounding medium, the important question 
arises as to whether the tunneling rate is enhanced or suppressed by the environment 
effect. There are already many works on this subject PP, and most past works deal 
with a system in equilibrium as a whole. The Euclidean technique such as the bounce 
solution [2j is often used in this context |SI-|H]- We would like to examine this problem 
by working in a new formalism of the real-time description of tunneling phenomena 
We find it more illuminating to use a real-time formalism instead of the 
Euclidean method much employed in the literature. An advantage of the real-time 
formalism is the possibility of identifying an enhanced tunneling rate at a finite time 
of dynamical evolution of the tunneling ^11 . 

Our formalism is based on separation of a subsystem from thermal environment, 
and integrating out the environment degrees of freedom. This method is best suited 
to a (by itself) non-equilibrium system which is immersed in a larger thermal equi- 
librium state. Detailed properties of the environment and its interaction form with 
the subsystem are expected to be insensitive to the behavior of the subsystem in 
question. We use the standard model of environment consisting of an infinite num- 
ber of harmonic oscillators jS], In this picture dissipation seen in the behavior 
of the subsystem is due to our ignorance of the huge environment degrees of freedom. 
Although the subsystem itself is an open system having an interaction with environ- 
ment, the entire system including the environment obeys quantum mechanical laws 
in the closed system. 

We work out consequences of the real-time formalism for a simple subsystem in 
one dimension. The most interesting result that comes out in this study is the reso- 
nance enhanced tunneling Jl]. The enhancement occurs at a time scale of order l/r/, 
with 7] the friction coefficient. The enhanced tunneling probability goes with t] like 
T]^^-^ according to our numerical analysis in the present work. We offer an interpre- 
tation of the resonance enhancement due to an energy flow into the system from en- 
vironment. Furthermore, the enhanced tunneling rate is related to the phenomenon 
of the parametric resonance. We call this enhanced tunneling as environment-driven 
excited tunneling (EET). The popular Euclidean approach cannot deal with the fi- 
nite time behavior of tunneling, hence misses the resonance enhancement at finite 
times. 
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Our method, with a certain extension, should be apphcable to field theory models. 
In future we hope to investigate applications to cosmology, in which one may deal 
with the electroweak first order phase transition relevant to baryogenesis JSj- An 
important question here is the time scale of tunneling. If the usual Hubble time 
scale that gives rise to the potential (or more properly the free energy) change, 
hence regarded in the conventional picture as the termination of the first order phase 
transition, is replaced by a shorter tunneling time scale of order l/r/, the conventional 
picture of the first order phase transition may drastically be changed. For instance, 
bubbles of the true vacuum may be formed by the resonance enhanced tunneling prior 
to the potential change. Since this process is stochastic, it may enhance the out- 
of-equilibrium condition necessary for baryogenesis J3], The bubble formation 
via EET may also resurrect the old scenario of GUT phase transition which was 
once discarded due to the graceful exit problem [T7] . 

The rest of this paper is organized as follows. In Section II the model and the 
basic tool of analysis is explained. In Section III the semiclassical approximation 
is introduced and the result of the crucial simplification, the linear approximation 
of the environment variable, is worked out. This approximation is used while the 
particle in one dimensional potential is in the metastable region, to determine how 
the particle is excited in the potential well due to the environment interaction. In 
Section IV detailed numerical analysis in the metastable region is presented, taking 
a concrete example of the potential. It is shown that excitation to higher levels 
proceeds via the parametric resonance effect in the linear approximation. In Section 
V the semiclassical formula in preceding sections is combined with the quantum 
mechanical barrier penetration factor to derive the tunneling rate in thermal medium. 
Numerical analysis is then presented, to quantitatively demonstrate the effect of 
EET. In the final section application to cosmological baryogenesis is briefiy discussed. 
In Appendix we give the friction coefficient along with relevant formulas in the 
standard model of particle physics, which may be applicable to discussion of the 
electroweak baryogenesis in the standard model. 

II Model and basic tool 

We consider as the first step of our investigation the simplest, yet the most basic 
problem of this kind, one dimensional system described by a potential V{q). The 
potential V{q), as illustrated in Fig^ is assumed to have some local minimum at 
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Figure 1: A typical potential for tunneling in one dimensional system. Curvature param- 
eters, Cor at the bottom of the potential, ujr at the top of the potential, along with the 
barrier height Vh and the barrier width, characterize a global structure of the potential. 

g = with V(0) = 0, which is separated at the barrier top, q = qs {> ^) , from a 
global minimum, which we take in the present work to be at g = oo such that there is 
no reflection from a wall at the far right. This system is put in thermal medium which 
acts as an environment. It has infinitely many, continuously distributed harmonic 
oscillators given by their coordinates Q{uj). Its coupling to the tunneling system is 
given by a Hamiltonian [3], [T^ . 



The frequency dependent coupling strength is c{uj) and Uc is some threshold fre- 
quency. One may imagine a generalized case in which the subsystem variable q is 
the order parameter for the first order phase transition in cosmology, the homoge- 
neous Higgs field, and the environment oscillator Q{uj) is a collection of various forms 
of matter fields coupled to the Higgs field. 

The basic dynamical equation in our problem is then 



One may eliminate the environment variable Q{(jj) from the second equation above, 
to get the Langevin equation [TH], 







(3) 



d q dv ^ /■* , . , . ^ , 



(4) 
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The kernel function aj here is given by 

POO 

ai{t) = — / dLor{ijj)sm{uot) , (5) 



Pl3{Q,Q') = 



exp 



iQ' + Q") cosh{[3uj) - 2QQ' 



with 

= ■ (6) 

By definition the spectral function r{uj) is positive definite. The random force term 
Fqlt) in eq.(j3]) is linear in initial values of environment variables, Qi{uj) and Pi{uj) = 
Qi{uj) at t = 0; 

FqH) = - r doocioo) f g,(cu)cosM + P,M^^^^ . (7) 

JuJc \ UJ j 

By taking the thermal bath of temperature T = 1/ (3 given by the density matrix, 

\ 1/2 

^TT coth(/5cj/2)/ 

UJ 

2 sinh(/5cj) 

for each environment oscillator, the following correlation formula in thermal medium 
is obtained; 

/■°° I3uj 
{{Fq{t) ,FQ{s)}+)env = duj r (u) COS UJ {t - s) coth{—) = anir - s) . (9) 

The combination of the two real a's, oijiif) + %ai{t\ is the well-known real-time 
thermal Green's function summed over oscillators, being added with the 

weight (?{uj^. 

An often used simplification is the local, Ohmic approximation taking 

r{uj) = -L- , (10) 

TT 

with ujc = 0, which amounts to fK)\ 

ai{T) = 6uj^6{T)+7]6'{T). (11) 

This gives the local version of the Langevin equation, 

d'^q dV ^ n dq ^ , . 

+ _ + + (12) 

The parameter 5uj'^ is interpreted as a potential renormalization or a mass renormal- 
ization in the field theory analogy, since by changing the bare frequency parameter 
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to the renormalized uj^ the term 600"^ q is cancelled by a counter term in the poten- 
tial. On the other hand, 77 is the Ohmic friction coefficient. This local approximation 
breaks down both at early and at late times [21], but it is useful in other time regions. 

The value of the coefficient of quadratic term of the potential is shifted by en- 
vironment interaction. We discuss the relation between the renormalized curva- 
ture ujr (oJr) at the barrier top (the potential bottom) and the physical curva- 
ture ub (i^p) here. These physical curvatures ujb,^p are frequently used in the 
following discussion. We ffist illustrate this for the curvature at the barrier top. 
In the vicinity of the barrier top, the renormalized potential is given by V{q) ~ 
—uj\{q — qsY /'^ + Vh + counter terms {ojr > by definition). The physical curva- 
ture ujb (the tachyon mass in the field theory analogy, which is defined by a positive 
value) shifted by the environment is determined as a pole position of the Green's 
function (0|T[g(i(:)g(0)]|0) at the barrier top, thus as a solution for the isolated pole 
of the following equation; 



In many past works, the renormalized curvature ujr was used as a reference 
curvature to evaluate the environment effect. Due to the positivity of r{ijj), one 
has uj'^Q < uj\ in general. A consequence of this is that the potential barrier is 
further suppressed by using ojb instead of taking the renormalized curvature ojr as 
the reference point. It is instructive to give an example. In the Ohmic model of 
eq. (jl(J|) the explicit form of the physical curvature cj^ is given by 



The relation between the renormalized and the physical curvature at the local min- 
imum of the potential is obtained in the same way, to give in the Ohmic model 



In subsequent discussion we shall specify which parameter we use. 

We now turn to description of the macroscopic system determined by dynami- 
cal coordinates q,Q{uj). Any mixture of quantum pure states, namely a quantum 
statistical state, is described by the density matrix. 




(13) 




(14) 




(15) 



00 



Wn \n){n 



(16) 



n=0 



6 



where \n) is a linearly independent state vector for a pure quantum state and < 
w„ < 1. As a special case one may deal with a pure state taking w„ = for n 7^ 0. 
The density matrix that describes the entire quantum mechanical system obeys the 
equation of motion; 

^/^§^ = [i^tot,p], (17) 

where iftot is the total Hamiltonian of the entire system. For the time being, we 
explicitly write the Planck constant h. In the configuration space this density matrix 
is given by its matrix elements, 

p(g,g'; QM , Q'H) = (g , QHI p |g' , Q'H) . (18) 

Its Fourier transform with respect to relative coordinates, q — q' , Q{uj) — Q\uj) , is 
called the Wigner function and is denoted by fw] 

/oo r 
n dX{uj) eM-^^P - * / dujX{uj)P{uj)] 

X p(g + i/2,q- e/2 ; Q{to) + X{uj)/2 , Q{lo) - X{lo)/2) . (19) 
It is easy to show from the master equation p7|l that this quantity obeys ; 

dfw dfw 1 / / , «^ ^ \ T. / ^ ^ 1 ^ fon\ 

^ + J|) + ^^ak) + ) ■ 

Reduction to the subsystem by integrating out environment variables Q{uj) is an 
essential part for our physical understanding of tunneling in thermal medium. We 
define the reduced density matrix of the tunneling system by 

p(^)(g,g')= / dQ{u;)dQ'{uj) 6{Q{uj) - Q'{uj)) ■ p{q ,q' ; Q{uj) ,Q'{uj)) . (21) 



The fact that there is a delta function 6{Q{uj) — Q'{uj)) for the environment variable 
reflects our understanding that one does not observe the environment. A simi- 
lar reduction is made to the Wigner function, to get the reduced Wigner function 

fi^\q,P) = J dP{to) dQ{u;)fw{q,p;Q{u;),P{u;)). (22) 

The reduction to the subsystem here is in spirit opposite to that common in 
the effective field theory approach; integration over heavy degrees of freedom is 



7 



performed. What is integrated out in the present case is the majority of constituents 
which is taken invariant under frequent thermahzing interaction. 



Ill Basic idea and semiclassical Wigner function 

In considering our problem of tunneling in thermal medium we assume that effects 
of thermal environment are most important in the time period when the particle in 
one dimensional system remains in the metastable potential well. Thus, we anticipate 
that at tunneling the simple quantum mechanical formula may directly be used. We 
thus separate the problem into two parts; first to treat the potential region left to 
the barrier top, g < g^, by semiclassically integrating out environment variables and 
then by projecting the derived reduced density matrix onto the energy eigenstate of 
the subsystem, which makes it possible to convolute the usual barrier penetration 
factor used in quantum mechanics. In the following two sections we do this and 
derive a formula for the tunneling rate in thermal medium. 

In the semiclassical ^ — > limit we have in eq. ()2U|) 

The resulting equation for the Wigner function, 

dfw ^ _ dfw_ ^ dV_ dfw_ ^24) 
dt dq dq dp 

is identical to the classical Liouville equation familiar in classical statistical mechan- 
ics. It has an obvious, formal solution; 



fw{q,P,Q ,P) = J dqidpi J dQidPi {qi ,Pi,Qi, Pi) 

■S {q - gel) S {p - Pel) S {Q - Qd) S {P - Pel) , (25) 

where qci,Qd{^^) etc. are the solution of (0), © taken as the classical equation 
with the specified initial condition written in terms of g^ , Pi 

1 Qi ) Pi- 

We consider the circumstance under which the tunneling system is initially in a 
state uncorrelated to the rest of environment. Thus we take the form of the initial 
Wigner function of a factorized form, 

/«(g, , , P,) = fi^\g{q^,p^) X /«^q(Q, , P,) , (26) 
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to get the reduced Wigner function after the environment Q{uj) ,P{uj) integration, 

fw\Q^P'^'t) = J dqidpi fw,g{qi , Pi) K{q ,p,qi,pi;t), (27) 
K{q ,p,qi,pi;t) = J dQidPi fw,QiQi , Pi) Hi - Id) Hp ' Pa) ■ (28) 

The quantity qd is now a solution of the Langevin equation (j3)), regarded as a function 
of initial values; qc\{qi ,Pi ,Qi , Pi]t). The problem of great interest is how further 
one can practically simplify the kernel function K introduced here. 

In many situations one is interested in the tunneling probability when the envi- 
ronment temperature is low enough. At low temperatures of T <^ a typical frequency 
or curvature scale a;^ of the potential V, one has 



UJ.« 



Q2(^),V^'M =0[Vf]«y^. (29) 



In the electroweak and the GUT phase transition this corresponds to T <^ m//(IIiggs 
mass). Expansion of qd in terms of Qi{uj) , Pi{uj) and its truncation to linear terms 
is then justified. Thus, we use 

dXn 



Hq-Qci) = f — ^ exp Ag ( g - gc 
J Zn 



dXg 
27 



a, [q-qS^-j dio{Q,iuj)q;^\uj) + P,iuj)qiPiuj)} 

valid to the first order of Qi{uj) , Pi{uj). The zero-th order term is q^^ obeying 
d'q^' 



, (30) 



dt^ 



+ (?-) +2 f dsajit-s)qS\s)=0, (31) 
\dq J (0) Jo 

^cl 

while the coefficient functions of the first order terms are q'^^\uj) , q^P (ui) obey- 
ing equations later shown, eq.(|T7|). A similar expansion for 5(p — Pci) using 
p^^\ p^^\uj) , pIi\u) , also holds. 

An alternative justification of the expansion (jHUj) is to keep 0[c(ti;)] terms consis- 
tently in the exponent since both q'^p{uj) and q^p{uj) are of this order. Thus, both in 
the low temperature and in the weak coupling case our subsequent analysis follows. 

With one can proceed further. Gaussian integral for the variables Qi{uj), 
Pi{uj), first, and then for Xg , Xp can be done explicitly with eq.()30|). when one takes 
the thermal, hence Gaussian density matrix for the initial environment variable: 



fw,QiQi ^Pi) = - tanh ^ exp 

71 Z 



- tanh ^ I ujQ] + -j- 



(32) 
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The result of this Gaussian integral leads to an integral transform ^01 of the Wigner 
function, f^) (initial) (final), using the following form of the kernel function; 



K[q ,p,qi,pi;t) = — — — exp 



1/ (0) {o)x ^ ( Q - Qci 

- 7;{(1 - (lci\ P ~ Pel ) J (0) 
^ \P-Pci 



, (33) 



where the matrix elements of J'^j ^ = lij are given by 



1 Buj 1 

hi{t) = - / du coth^-\z{u;,t)\\ (34) 

1 du) 1 

/22(t) = - / duj cotYi- |z(cj,t)p, (35) 



2, J 2, u! 

h2{t) = l r dLucoth^-?ftiziu,t)z*iiu,t)) = ^-^, (36) 
Z Jujc Z UJ Z 

where z{uj ,t) = q^.^^ {uj ,t) +iuj q^^"^ {uJ ,t) , and z{uj ,t) = p^f ■* {u ,t) + iuo p^f ^ (c^ , t) • 
Quantities that appear in the integral transform are determined by solving dif- 
ferential equations; the homogeneous Langevin equation for q^^ , eq.(|^ and an 
inhomogeneous linear equation for z{uj ,t) and z{uo ,t), 

d'^ziu f) / d'^V\ /■* 

a zyiv, ) ^lav^ ;2(cu , t) + 2 / ds ai{t - s)z{uj , s) = - c(cu)e*"* . (37) 



dt'^ \ dq^ / (0) JO 
A similar equation as for z{uj ^t) holds for z{(jJ ,t) . The initial condition is 

qf{t = {]) = q,, p'^\t = 0)=pi, zico,t = 0) = 0, ziu,t = 0) = 0. (38) 

The physical picture underlying the formula for the integral transform, eq. (j27|) 
along with (jSS)), should be evident; the probability at a phase space point {q,p) 
is dominated by the semiclassical trajectory q^,^^ (environment effect of dissipation 
being included in its determination by the third term of eq.(j2^) reaching (g ,p) from 
an initial point (gj ,pi) whose contribution is weighed by the probability fl^ initially 
given. The contributing trajectory is broadened by the environment interaction 
with the width factor ^Jhj- The quantity In given by (jH^ . for instance, is equal to 
{q — ql,i^y ; an environment driven fluctuation under the stochastic force FQ{t). 

We note that in the exactly solvable model of inverted harmonic oscillator po- 
tential the identical form of the integral transform eq. ()33|) was derived [Oj without 
resort to the semiclassical approximation. Explicit form of the classical and the fluc- 
tuation functions {q^i^ and z{uj ,t)) is given there, which is valid beyond the Ohmic 
approximation . 
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It would be instructive to write solutions in the harmonic approximation taking 
the frequency at the local potential minimum with LUp = \J uj\ — ^ , eq. ()15|) . the 
potential is given by V{q) = ^uj^q^. Then the zero-th order term q^{t) and z{ijj,t) 
are obtained analytically, and given by 

qf =u{t)qi + v{t)pi, (39) 

u{t) = icosujpt + ^smUpt)e-^'/\ vit) = !!^i^e-^*/^ (40) 
2,ujp ujp 

z{uj , t) = „ ""^^^ . ( e'^' - u{t) - luvit) ) . (41) 

We used the local Ohmic approximation such as 

ft 
/o 



dsaiit — s)z(oJ ,s) — > r] z{lj ,t) . (42) 



The function z{uj , t) has a Breit-Wigner form for a small rj. In particular, the infinite 
time {t oo) limit of the fluctuation has 

/n(oo) = du; coth - • (43) 

The friction t] appears as the resonance width in this formula. The real part of 
the pole position in the frequency integral such as PHj). Up = ^ (pole position), is 

The harmonic approximation is excellent at late times, but at early and inter- 
mediate times nonlinear resonance effects are important as discussed in [11.. The 
nonlinear effect will be taken into account in discussion of the finite time behavior 
of the tunneling rate. 

A convenient form of the fluctuation formula is derived by separating the zero- 
temperature result along with coth/?L<j/2 = 1 + 2/(e^'^ — 1). We denote the tem- 

(T) (0) 

perature dependent part of the fluctuation by I^j . The zero-temperature part I^j 
must be renormalized. After performing renormalization the zero-temperature part 
is approximated by result of the harmonic approximation. The renormalized fluctu- 
ation is then well approximated using the physical curvature, ujp, eq. pSj) . We shall 
use the same notation lij for the renormalized fluctuation in the rest of our paper. 



IV Numerical analysis in the metastable region; relevance of para- 
metric resonance 
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Figure 2: The tunneling potential used in numerical calculation. 

We point out a relation of the semiclassical Wigner function in the preceding 
section to the parametric resonance, as briefly mentioned in [llj. Presence of an- 
harmonic terms is inevitable in any realistic tunneling potential. This introduces 
a non-trivial oscillating term in the coefficient function (^r) (o) of the fluctuation 
equation, (jHTj) . thus giving a differential equation akin to the Mathieu equation. 
Depending on whether or not the relevant parameters fall into the instability band 
of this Mathiew-type equation, the fluctuation may increase indefinitely in the zero 
friction limit fl^. What happens is more subtle; the parameters fall right on the 
boundary of stability and instability bands. We shall show this by taking an example 
of the tunneling potential. 

The potential we take (illustrated in Fig|2)) is given by 

V{x) = -^x^ + guj\x^, {x < Xb) 

= Vh, {xb<x<Xf), =--1/(00), {xf<x). (44) 

The two parameters Vh and x^ are related since we took a particular quartic potential 
in the metastable well. 

The fluctuation z{uj,t) in the Ohmic approximation follows the equation; 

z{u;, t) + V"{qf) z{uj, t) + r] z{u;, t) = - c{lj) e'^' . (45) 

The potential derivative here is 

V"{cif)=^l + l2gCol{qff. (46) 

Taking the zero friction limit 77 = gives the following related homogeneous equation, 

i{uj,t) + [h + 2e{q^^)M^.t) = Q, (47) 
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Figure 3: Band structure for relevant fluctuation. Shaded areas correspond to the in- 
stability band. 1st band is enlarged separately (the left figure), where the crossed point 
corresponds to our model relation, ea.()48|). 

where 

h = ul, e = QgCul. (48) 
The equation for q^^^ is as follows: 

+ V'iqS^) = (49) 

This gives a periodic solution for A notable feature of the homogeneous equation 
in the rj = limit (|47p is that q^^^ appearing in the coefficient function is periodic, 
hence it gives rise to solutions either of the Bloch wave type or of the parametric 
resonance type [22] • Despite of the definite relation, eq. ljlHj) . we regard {h,9) as 
free parameters to throughly investigate the stable-unstable band structures, be- 
cause that would more definitely clarify relation between the band structure and the 
relevant parameter region of this model. 

For simphcity we analyzed the problem assuming 

g = OM, Vh = hup, Up = ^Col-7]^/A . (50) 

In a relevant region of 

Q<e <2Cj\, 0<h<5ul, (51) 

which includes the model region ()48|) . we numerically checked the time evolution 
of the homogeneous solution of eq. ()47p . We observed the power-law increase of the 
fiuctuation. Indefinite increase of the fiuctuation z{uj,t) as time increase implies 
that the assumed parameter lies either within the unstable band or on the boundary 
between the bands. The linear rise z{uj, t) oc t, as indicated by our numerical analysis 
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Figure 4: Time evolution of the function y{t), a solution of ea. 1)54(1 . The amplitude of the 
function y{t) with the friction rj = O.OOSwp is depicted. 

suggests that the parameter is on the boundary. The result of the band structure 
is given in FigOl and the relevant parameter eq. ()48p (corresponding to the crossed 
point) is just on the boundary between stable and 1st unstable bands. This is 
consistent with that the linear rise, as is expected in the case of the exact Mathiew 
function on the boundary 

For this numerical calculation we had to choose an initial condition for q''^\ We 
took several conditions for investigating band structures, and we always obtained 
the same result of the linear power-law increase. 

We next consider effects of the friction rj. The inclusion of the friction gives a 
modification of eq. ()47|) to 

Zhit) + [h + 29 iqTf]zh{t) + 7]Zh{t) = , (52) 

The effect of the friction term riZh{t) is studied by changing the variable Zh{t) to yit), 

y{t) = z,{t)e^'/' . (53) 

Then the equation for y{t) becomes 

y+[(h-r^yA) + 29{q^^f]y = 0. (54) 

On the other hand, the equation for q*^^ is 

qS^ + V'iqS^)+mci{0) = 0- (55) 

When the effect of friction for the behavior of the zero-th order q^'' is small, the 
equation for y is nearly of the parametric resonance type. In this case then, the 
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Figure 5: Examples of u) dependence of z{ijj,t). The quantity ujp\z{uj ,t)\'^ + \z{u} ,t)\^ / u)p 
both a) at early and b) late times are shown. The unit of time is ujp"^- 

behavior of z{uj,t) is rather simple; it is a product of the linearly rising function 
{yit) ~ t) and the exponentially decreasing function (e"''*/^). The rate of the expo- 
nential decrease is 77/2, hence this product function has a maximum at a time around 
2/7]. This seems essentially what we observe in numerical computation for z{uj,t). 

With the presence of the friction the function y approaches the behavior of simple 
harmonic oscillator as time increases since g^^^ by friction. This means, as seen 
in Fig^ that the linear increase of y is saturated at the time of order 2///. The net 
effect of the friction is simply the minor change of a coefficient of A/rj, where A is 
close to 1.7. 

The effect of the presence of the inhomogeneous term in eq. ()45|) is as follows. 
This term acts as an external force, and gives a resonance effect when the frequency 
uj is close to the inherent frequency of the system. We illustrate this point. In 
FiglSlciJ-dependence of ujp\z{(jj,t)\'^ + \z{(jj,t)\'^ /ujp at a particular time is shown, and 
one observes, besides the simple resonance, higher modes due to the nonlinear effect. 
Without the strong resonance effect like this one, the fluctuation z{uj , t) is very small, 
since we have the initial condition, ^(07,0) = = i(co',0). Thus, the environment 
acts here as a force of resonant kick from the zero amplitude, which then makes 
possible the parametric resonance effect to work, typically very important for large 
amplitude oscillation. In FigEK), the growth of 1st {u ~ Up) and 2nd {u ~ ?>ujp) 
mode is shown. Other higher modes (3rd, 4th, ■ ■ ■) are also amplified in the same 
manner as in the 2nd mode case at t ~ 360ti;~^. On the other hand, these higher 
modes disappear at late times, and only the first mode remains as seen in FiglHb). 
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Figure 6: Time evolution of a) the fluctuation u>pl\i + /^p ^^'^ b) the subsystem 
energy {Hsys.) - (T = part). A temperature T = 0.5ujp is taken with four cases of the 
friction ij (in the unit of ujp). 

Thus, the growth of higher modes is a finite time effect. 

We now discuss the behavior of the frequency integrated fluctuation I{t) in 
eqs. fj34p - ()36|) . The increase of the fluctuation due to the parametric resonance gives 
an interesting time evolution of I. In FiglHK); the time evolution of the quantity 
cijp/^P + I^P l^p is shown for a few values of the friction rj. This function depends 
on the initial value, Qi ,Pi, for which we took a typical value in the ground state at 
the potential minimum. 

What is the physical picture behind the resonance enhancement ? We would 
like to present a suggestive interpretation. We first derive a relation between the 
fluctuation here and the subsystem energy defined by 

(Hsys.) = (y + y^ci)- (56) 

By making the approximation on g^i and Pc\, namely taking up to linear terms of the 
initial environment value as used in eq.ipUll. this equation is reduced to 

(Hsys.) = (T = part) 

Thus, the temperature dependent part of the subsystem energy is given by the 
fluctuation averaged over initial values. In short, the subsystem energy is driven by 
the fluctuation. If initial values for (gci(O)); (Pci(O)) are small (in the Up = 1 
unit), the zero-temperature part (T = 0) part is close to 0.5u!p (zero-point fluctuation 
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of the system). In FiglHb) the quantity Hgys— (T = part) is shown, taking as the 
initial state the ground state of the harmonic oscillator at the potential bottom. The 
ground state we use is characterized by 



The time at which the maximum of H^ys occurs is 1-7 /rj, which is the same as for 
that of the fluctuation, thus indicating a close correlation between the fluctuation 
and our subsystem energy {Hgys.)- We shall have more to say on an important role 
of (Hsys.) towards the end of the next section. 

Implication to the tunneling rate seems obvious; the parametric resonance we 
observe here enhances excitation to higher energy states at the potential region 
left to the barrier top in the quantum mechanical terminology. This would then 
enhance the tunneling rate from the metastable state of the potential well. We shall 
demonstrate this in the following section. 

V Tunneling rate formula 

Energy eigenstates are special among pure quantum states, since they time-evolve 
only with phase factors. Thus, if one takes for the projection in eq. (jl6p the 

energy eigenstate \n) of the total Hamiltonian, the density matrix does not change 
with time. On the other hand, if one takes the energy eigenstate of the subsystem 
Hamiltonian, then the density matrix changes solely due to the environment interac- 
tion. It is thus best to use a pure eigenstate of, or its superposition of, the subsystem 
Hamiltonian, when one wishes to determine how the tunneling rate is modified in 
thermal medium. 

We thus perform projection onto the energy eigenstate of the tunneling system. 
The projected probability onto energy eigenstate of the subsystem Hamiltonian is in 
general given by 



where ^e{<1) is the wave function of eigen energy E. 

Our basic idea is to limit application of this formula to the left region of the 
metastable state and to use the semiclassical formula for the Wigner function 



fwili^Pi) = -exp -ujpQi 




(58) 



p;t), (59) 
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fwiQyP ',t) derived in the preceding section. Explicit p integration in gives 

1 



exp 



27r 

1 ( q + q' 



2/11 



-qf{t) 



det/ 



i\2 



(o)/,N , h2 ,q + q' (0) 



+ ^{q'-q) W + y^^. ^ 



2/11 
-q'^{t)) 



q 



(60) 



The tunnehng probabihty in medium is then obtained by multiplying the well-known 
formula of the barrier penetration factor. 

Unless the energy eigenstate near the barrier top is important for computation 
of the tunneling rate, one may approximate the wave function in the left region to 
q = qs by that of the harmonic oscillator at the potential bottom, g = 0. The wave 
function fEiq) is thus approximated by using the Hermite polynomial Hn as 



^E{q) 



9V2 



(61) 



V7r2"n!' 

where Up is the shifted frequency at the local potential minimum, eq. ()15|l . 

To compute the projected probability VVn(t), we introduce the generating func- 
tion; 



w(t,z) = £ w„(t);2r 



(62) 



n=0 



The following formula of the Hermite polynomial, the Mehler's formula 



E 

n=0 



[z/2y 



2xyz — (x^ -I- y'^)z'^ 
1-^2 



(63) 



is useful. From eqs.flUnj). (jUTj) . (jU^ . one has Gaussian integrals for q and q', resulting 
in 



Wit,z) 



A 



27r 



fwiqi^Pi) 



1 



det A 



- exp 



1 - z 



-Aqf.vt^)A 



2 det A 



(0) 

(0) 
Pel 



(l-^)(J22-f)+a;p 



:i - z)i 
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-(1-.)/ 



12 



;i--)an-2^) + :^ 



(64) 
(65) 



Quantities Jjj are to be computed at finite time t. The important constraint of 
the probability conservation Y^^=q H'„(t) = 1 is automatically satisfied due to 
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Figure 7: Time evolution of the projected probability for the ground, the second excited, 
and the forth excited state. The case of T = O.biVp and r] = O.OOSwp is depicted. For 
comparison the result in the harmonic approximation is shown. 

W{t ,1) = 1. For precise computation, one expands this generating function to 
powers of z, and identifies the projected probability yV„(t) by (jU^ . 
We next clarify the relation, 

oo 

EnWn{t) = {Hsys.) (66) 

n=0 

where {Hgys.) is calculated by ^7\. This relation was numerically checked to be 
correct, as might be expected. We have no analytic proof of this relation, although 
there is no doubt on the correctness of this relation. 

All numerical calculations were performed by using the concrete example of the 
tunneling potential eq.(jl3)) throughout this section. We assume the subsystem is 
in the ground state at the initial moment, and use eq. ()58|) as the initial Wigner 
function. Some results were compared with that of the Harmonic approximation 
discussed in the last of the section III. The effect of the non-linear term is clarified 
by this comparison, because the term is neglected in the Harmonic approximation. 

In Figl?! a few examples of the projected probability are shown, assuming the 
initial ground state; >Vn(0) = 6no- It is clearly observed that excitation to higher 
energy levels is correlated to evacuation of the ground state. The infinite time limit 
of yVn{oo) coincides with the result of the harmonic approximation at the potential 
bottom, 

W„(oo) = (1 - e-^'^)e-"^- + 0{ri) , (67) 
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Figure 8: a) Time evolution of th tunneling rate T(t). A temperature T = O.Swp is taken 
with four cases of the friction r]. b) Ratio of the maximal tunneling rate at t ~ 0{2/rj) to 
the r(oo) (a traditional one) is depicted as a contour map in the rj — T plane. 

as is expected. We note that the subsystem thermalizes at the infinite time de- 
spite that it was initially in the ground state of a harmonic oscillator, hence a non- 
equilibrium state. 

The tunneling rate r(t) per unit time is roughly the tunneling probability , 
since ujp/27c = 1 /period for the periodic motion in the left potential well. Thus, 

r(t) = /E>v„(t)im)p, (68) 

2^ n=0 

e{Vh -E)+ e{E - Vh) . (69) 

In computing the tunneling rate we use a modified potential V which is obtained 
from V by replacing the curvature parameter ojr hy ujb = \J^lk + ff — ri/2. The 
reason for this is that one should compare the tunneling rate with and without 
the environment effect by using effectively the same values of parameters for the two 
cases. This modification explains the bulk of the environment effect, the environment 
suppression factor proposed by |3j, as shown in p. In equations above qi{E) are 
turning points separating the subbarrier region. 

In FiglH^) the tunneling rate r(t) is shown for a few values of the friction. The 
tunneling rate is observed to become maximal at the time t ~ ^-T/v- salient 
feature, and the most important result of the present work, is that the tunneling rate 
is enhanced around the time of order 2///, which is caused by non-linear resonance 
effects. The enhanced tunneling rate at t ~ 0(2/?7) is clearly related to the power- 



|T(E)|^^exp -2 / dxJ2(V(x) - E) 
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law rise of the fluctuations la, as observed in our previous paper [TT] and as already 
explained. The maximal enhancement factor for the tunneling rate is plotted in 
FiglHb), relative to the harmonic approximation. A large enhancement is obtained 
for small f] and small T values. In the harmonic approximation an almost monotonic 
behavior of time evolution is evidently observed, thus the rate being given by r(oo). 
We may summarize what happens, in the following way. The non-linear resonance 
inherent to the tunneling potential excites higher modes that have larger tunneling 
probabilities. This effect is however terminated by the friction caused by environment 
interaction, which is the origin of the final decrease of the tunneling rate at large 
times. 

We finally study at a quantitative level correlation between the tunneling rate 
r{t) and the tunneling system energy Hsys{t). In order to evaluate the correlation 
quantitatively, we define the correlation coefficient r used in statistics; 



{V){Hsys) ' Time Jo 

The correlation coefficient r was numerically calculated and we obtained results that 
can be summarized by r = 0.99 ~ 1 in the parameter region of temperature and 
friction, 0.2 < T/ujp < 1.0 and 0.001 < rj/uop < 0.02. This proves the correlation 
between the two quantities. 

We next turn to the discussion of the survival probability, the probability that 
the metastable state remains in the same initial state. This quantity P{t) is defined 

by 

p(t) = -r(t)P(t), P(o) = i. (71) 

This gives 



P{t) = exp 







dt' T{t') 



(72) 



In FiglHti') the survival probability is depicted along with the result of the harmonic 
approximation of the potential bottom. The effect of EET is unquestionable. 
A traditional definition of the lifetime ttrad which many authors use is 

ttrad = r^^(oo) , (73) 

while in the case of time dependent tunneling one may effectively use ^lt given by 
the e-holding time, 

P(tLT) = . (74) 
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Figure 9: a) Time evolution of the survival probability P{t). The case of T = O.bUp 
and rj = O.OOScUp is depicted. For comparison the result in the harmonic approxima- 
tion is shown, b) Ratio of the lifetime as defined in eq. ()74|) to the traditional one is 
depicted as a contour map in the rj — T plane. 

In FiglHb) a contour plot of tLT/^trad. is shown in the parameter ij] ,T) space. For 
small values of r] and T a substantial reduction of the lifetime ^lt is observed, which 
means the early termination of the tunneling. 

We hope that we have convinced the reader of the power and the usefulness 
of EET, although our idea is vindicated by detailed numerical analysis which is 
necessarily limited to special cases. 

VI Discussion 

There may be many applications of environment- driven excited tunneling (EET) 
presented in this paper. We shall only mention two possible applications to cosmol- 
ogy, an enhancement mechanism of electroweak baryogenesis and a possible resur- 
rection of the old inflationary scenario of the GUT phase transition . Furthermore 
in the present discussion we shall focus on one aspect of the baryogenesis condition; 
the out-of-equilibrium condition ^1], Needless to say, there are many problems 
to be solved in actual application of these ideas to realistic models. 

Electroweak baryogenesis is expected to occur if the electroweak phase transition 
is of a strong first order. The first order phase transition usually proceeds via the 
bubble formation. Inside the bubble the true electroweak broken phase is realized in 
which the baryon number is effectively conserved. On the other hand, outside the 
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bubble the universe remains in the high temperature symmetric phase in which the 
baryon number is not conserved via sphaleron related pro cesses [22]. This mismatch 
inside and outside the bubble gives rise to a possibility of electroweak baryogenesis, 
although details of a promising mechanism are yet to be worked out. Thus, the 
strong first order of the phase transition has been considered to be a requisite for 
the out-of equilibrium condition necessary for the baryogenesis. 

In the ordinary circumstance, namely in the standard model, the strong first order 
phase transition requires a small Higgs mass, which is inconsistent with observation. 
Along with the small CP violation effect of Kobayashi-Maskawa theory, this excludes 
)26j the possibility j2ZI of the standard model as a viable model for baryogenesis. 
Thus, a more complicated model such as two Higgs doublet model needs to be invoked 
for the electroweak baryogenesis. Regardless of whether an extended Higgs model 
is needed, we would like to point out that the environment-driven excited tunneling 
(EET) found in this paper has a chance of enhancing the out-of equilibrium condition 
necessary for the baryogenesis. 

Although not much discussed in the literature, there is an important factor to 
consider for the out-of equilibrium condition. This is the factor Ttn, where F is the 
nucleation rate of the bubble and tfj = 1/H is the Hubble time. In the usual scenario 
of the first order phase transition the effectiveness of electroweak baryogenesis is in 
proportion to this quantity. The nucleation rate F in the usual estimate without EET 
is exponentially suppressed by the potential barrier, and the factor Ttn is much less 
than 1, which means that the first order phase transition is never completed by a 
merger of nucleated bubbles. Under this circumstance it is expected that the phase 
transition is terminated by proceeding to disappearance of the local minimum of 
the symmetric phase so that the condition for the first order phase transition is not 
met. The effectiveness of the out-of equilibrium condition is then reduced by a factor 
proportional to the quantity Ttn- 

On the other hand, with the aid of EET, the factor Ft may become of order unity 
at some finite time, and the phase transition may be completed much prior to the 
Hubble time, as discussed in previous sections. One should check in comparison to 
the Hubble time the new time scale 2/ri which is roughly the time at the maximum 
for the tunneling rate. The factor Ft thus may become of order unity at t ~ 2///. In 
this case there is no reduction factor like Ttn- 

We shall check numerical relation between two time scales; and 2/1]. The size 
of the friction t] in the standard model is of the order of the Higgs mass tuh times 
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some power of coupling factors, while the Hubble rate is given by 



H = 1/tH 



47^2 7-2 



= 106.75. 



(75) 



45 rupi ' 



Calculation, as summarized in Appendix, indicates that 




r, ^ ^a;,tanh(-^), cul = V"{0) ^ 0.2 x (T^ - - 0.04.;^) , (76) 



where yt is the top Yukawa coupling, and V{(j)) is the effective potential of the 
standard model at finite temperature, with v the vacuum expectation value of the 
Higgs field. Thus, the Hubble rate is much smaller than rj by the inverse Planck 
factor since — <^ yf. 



One may expect that even in the case of weak first order phase transition the 
quantum tunneling enhanced by EET may give rise to a sufficient condition for the 
out-of equilibrium in the early phase of the phase transition. Much remains to be 
seen with this regard. 

Another interesting effect of EET concerns resurrection of the old inflationary 
scenario of the GUT phase transition JHl- The old inflationary scenario was once 
rejected because the bubble formation proceeds too slowly, hence leaving a very in- 
homogeneous universe ^7j. The essential reason of this problem, called graceful exit 
problem, is that the exponentially suppressed tunneling rate is too small such that 
bubbles of true vacuum do not merge sufficiently to give rise to a homogeneous uni- 
verse which the inflationary scenario intends to derive. Our mechanism of enhanced 
tunneling may help this situation. The problem must however be addressed quanti- 
tatively; one has to obtain an early termination of the first order phase transition to 
promote merger of bubbles, yet has to obtain an sufficient inflation to solve intended 
problems. We leave this quantitative study to future. 
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Appendix 



A-1 Effective potential and critical temperature 

We shall first give the finite temperature effective potential of the standard 
model at 1-loop order, assuming that the Higgs mass is much smaller than other 
particle masses of the standard model such as the top quark, the weak boson, thus 
rriH << rrit , mz , mw- This is taken as a necessary condition for the first order phase 
transition. This condition however is not met in reality. Nevertheless, we take this 
case as an important illustration for more complicated cases. 

The zero temperature part of the effective potential [2H] is as follows; 



where v = 2A6GeV is the vacuum expectation value (VEV) of the Higgs field, and 
rria is the mass of a particle species a. The finite temperature effect is then given by 
the free energy [221, 



The Higgs contribution has been neglected, due to the assumption of tuh <^ masses 
of other particles. 

According to we may make the high temperature expansion of the effective 
potential with T <^ i;, to get 




(77) 




(78) 




(79) 



where the dimensionless coefficients are given by 




2m^ + m1 + 2mj 



(80) 



— Ami 
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Here as = 21og47r — 27 ~ 3.91, ap = 21og7r — 27 ~ 1.14. The critical temperature 
is computed from these, to give 

^0 _ ^0 



l-Ey{X{T)D) l-{8GeV/mHy 

At the value of Tc two local minima coincide, while Tq indicates where the symmetric 
phase of = disappears as the local minimum. It is thus important to note that 
the condition > To for the first order phase transition is met. 

A-2 Friction coefRcient in the standard model 

The friction coefficient t] is computed using the method of [3n],[3I]. According to 
j31j . it is possible to identify the two-body bilinear operator of standard model fields 
as the environment variable Q{uj). The most dominant contribution to the friction 
7] comes from the top loop diagram where the top- Yukawa coupling is of the form, 
yt(f)it/\/2. Thus, the Langevin equation of motion for the zero-mode Higgs field 0o 
becomes 

00 + V\(f)o) + 2 /* dsQa{t - s)Ms) = F{t) , (82) 
Jo 

({F(t),F(s)})e„v = 5Ra(t-s), (83) 
aixo) = ^J rf=^x(tt(x)tt(0))env = / dcu^^^e'^^^^ . (84) 

The real-time Green's function {tt{x)tt{0))cnv can be obtained by the analytic con- 
tinuation of the imaginary-time Green's function and it is easy to calculate the 
imaginary-time Green's function by using the Matsubara formalism. The necessary 
spectral weight r^iuj) is given by 

= 1/(0^ - - 2"-^ ' 

where p = |p| is the momentum magnitude and Up = (e^^ + 1)^^, assuming that the 
top quark is massless in the symmetric phase around T ^ Tq. 

In the vicinity of the false vacuum, namely ~ 0, the behavior of the zero- 
mode Higgs field exhibits the exponentially dumped oscillation. The solution of the 
homogeneous zero-mode Higgs field equation (j82|) is approximately obtained by 

rniuj) sin(co't) 



Mt) ^ 2 duj 



(a;2 - a;2) + ■n'^rjj{u) 

'^^T~2 — 2 2/ \ = sm c^pt exp 1) , 86 
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where Up is the curvature of the false vacuum {Up = V^"(0)), and the friction coeffi- 
cient in this case is computed from eq. ()85|) . The result is 



-rK) = — ^,tanh^ • (87) 
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